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Abstract 

We present a new measurement of the decay B~ t~V-j- with a semilcptonic B tagging method, 
using a data sample containing 657 x 10^ BB pairs collected at the T(4S') resonance with the 
Belle detector at the KEKB asymmetric e+e" collider. A sample of BB pairs are tagged by 
reconstructing one B meson decaying semileptonically. We detect the B~ t^V^- candidate in 
the recoil. We obtain a signal with a significance of 3.8 standard deviations including systematics, 
and measure the branching fraction to be B{B^ t~Vt-) = (1.65lQ'3y(stat)^Q3^(syst)) x 10^^. 
This result confirms the evidence for B~ t~Vt obtained in the previous Belle measurement with 
a hadronic B tagging method. 

PACS numbers: 13.20.-v, 13.25.Hw 
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The purely leptonic decay B~ t~V [JJ is of particular interest since it provides a direct 
measurement of the product of the Cabibbo-Kobayashi-Maskawa(CKM) matrix element 
Vub j2| and the B meson decay constant fs- In the Standard Model(SM), the branching 
fraction of the decay B~ r^V is given by 

B(B- r-V)sM = ( 1 - 4V/IIK..|V«, (1) 



Stt V m 



where Gi;' is the Fermi coupling constant, m,- and are the r lepton and B meson masses, 
and tb is the B~ lifetime. The dependence on the lepton mass arises from helicity conser- 
vation, which suppresses the muon and electron channels. Observation of B~ —>■ t~V could 
provide a direct measurement of fs- Physics beyond the SM, such as supersymmetry or 
two-Higgs doublet models, could suppress or enhance B{B^ t^v) to levels several times 
as large as the SM expectation through the introduction of a charged Higgs boson [3|,0]- The 
expected SM branching fraction from other experimental constraints is (0.93to;i2) x 10~^ jsf. 

The previous Belle measurement Q reported the first evidence of B~ — >■ t~V decay with 
a significance of 3.5 standard deviations (a), and measured the branching fraction to be 
B^B" T^Vr) = (1.79^o;|9(stat)lo;5i(syst)) x 10""^, using a full reconstruction tagging 
method. The BaBar Collaboration has reported a search for B^ —>■ t~V decay with hadronic 
tagging [3] and semileptonic tagging [8|] using 383 x 10^ BB pairs. They report a 2.6 a 
excess, combining the two measurements. No statistically significant enhancement relative 
to the SM expectation has been observed in previous experimental studies. To establish the 
B^ T~V signal, test consistency with the SM and search for a charged Higgs boson effect, 
we need more statistics. In this paper, we present a new measurement of B~ — >■ T~Vr from 
the Belle experiment with a semileptonic tagging method. 

We use a 605 fb~^ data sample containing 657 x 10^ BB pairs collected with the Belle 
detector at the KEKB asymmetric energy e^e~ (3.5 on 8 GeV) collider operating at the 
T(4S') resonance (a/s = 10.58 GeV) [9]. We also use a data sample of 68 fb~^ taken 
with a center of mass energy 60 MeV below the nominal T(4S') (off-resonance) mass for 
a background study. The Belle detector is a large-solid-angle magnetic spectrometer that 
consists of a silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an array 
of aerogel threshold Cherenkov counters (ACC), a barrel-like arrangement of time-of-fiight 
scintillation counters (TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) 
crystals located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. 
An iron fiux-return located outside of the coil is instrumented to detect K\ mesons and to 
identify muons (KLM). Two inner detector configurations were used. A 2.0 cm beampipe 
and a 3-layer silicon vertex detector was used for the first sample of 152 xlO^-B-B pairs, 
while a 1.5 cm beampipe, a 4-layer silicon detector and a small-cell inner drift chamber were 
used to record the remaining 505 xlO^-BS pairs The detector is described in detail 

elsewhere ll[ . 



We use a detailed Monte Carlo (MC) simulation based on GEANT [1^ to determine the 
signal selection efficiency and study the background. In order to reproduce effects of beam 
background, data taken with random triggers for each run period are overlaid on simulated 
events. The B^ t^v^ signal decay is generated by the EvtGen package [l3|. To model 
the background from e^e~ — > BB and continuum qq (g = m, c?, s, c) production processes, 
large BB and qq MC samples corresponding to about four times the data sample are used. 
We also use MC samples for various rare B decay processes such as charmless hadronic, 
radiative, electroweak decays and b ^ u semileptonic decays. We also use a MC sample for 
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r+r events. 



The strategy adopted for this analysis is same as in the previous measurements. We 
reconstruct one of the B mesons decaying semileptonically (referred to hereafter as -Btag) 
and compare the properties of the remaining particle (s) in the event (-Bsig) to those expected 
for signal and background. In order to avoid experimental bias, the signal region in data is 
not examined until the event selection criteria are finalized. 

We reconstruct the _Btag in B~ D*^i~i7 and B~ — >• D^i~i7 decays. For reconstruc- 
tion, we use D*^ —>■ D^tt^ and D'^'j decays. mesons are reconstructed in K~7t~^, K'tt^tt^ 
and fr~7r+7r^7r+. For -Bsig, we use r~ decays to only one charged particle and neutrinos: 
l^ViVr and ir^Ur- We require that no charged particle or 7r° remain in the event 

after removing the particles from the -Btag and -Bsig candidates. 

Charged particles are selected from well measured tracks reconstructed with the CDC 
and SVD originating from the interaction point. Electron candidates are identified based on 
a likelihood calculated using dE/dx in CDC, the response of ACC, ECL shower shape and 
the ratio of the ECL energy deposit and the track momentum. Muon candidates are selected 
using KLM hits associated to a charged track. Both muons and electrons are selected with 
efficiency greater than 90% in the momentum region above 1.2 GeV/c, and misidentification 
rates of less than 0.2%(1.5%) for electrons (muons). After selecting leptons, we distinguish 
charged kaons from pions based on a kaon likelihood derived from the TOF, ACC, and dE/dx 
measurements in the CDC. The typical kaon identification efficiency is more than 85% and 
the probability of misidentifying pions as kaons is about 8%. Photons are identified as 
isolated ECL clusters that are not matched to any charged track. 7r° candidates are selected 
from pairs of photons with invariant mass between 0.118 and 0.150 GeV/c^. Photons from 
7r° candidates used in meson reconstruction and photons from D*^ -D°7 decays are 
required to pass the following energy requirements: 50 MeV for the barrel, 100 MeV for 
the forward endcap and 150 MeV for the backward endcap. For low momentum vr^'s from 
-D*° —>■ D^Ti^ decay, we require the photon energy to be greater than 30 MeV. 

-D° meson candidates are selected from combinations of charged kaon, pion and 7r° can- 
didates. We require the invariant mass of candidates to be in the range 1.851 < Mjjo < 
1.879 GeV/c^ for B)0 ^ K-7i+ and K-n+n-7r+ decays, and 1.829 < M^o < 1.901 GeV/c^ 
for —>■ K^TT^TT^ decay. D*^ candidates are selected by combining the candidates 
with low momentum tc^ candidates and photons. For D*^ candidates, we require the mass 
difference - M^o to be in the range 0.1389 < - M^o < 0.1455 GeV/c^ and 
0.123 < - Mdo < 0.165 GeV/c^ for B)*° ^ B)°7r° and D*° D^j decays, respec- 
tively. The regions correspond to 3a from the nominal -D° mass or from the nominal 
mass difference. We select -Btag candidates using the lepton momentum P/ and the co- 
sine of the angle between the direction of the -Btag momentum and the direction of the 
momentum sum of the D^*)" and the lepton cos^5_£,(*)£. This angle is calculated using 
cos9b_d(*)(, = (2Ebeam^D{*)£ " "^B " ^]^(*) () / {'^Pb ' PdMi), whcrc Ejjc,)^, Pd{*)i and Mj^(,)^ 
are the energy sum, momentum sum and invariant mass of the and lepton. All pa- 



the momentum of B meson in the cms calculated with Pb = y E^^^^ — m^. For the signal 

side track, we require the momentum P*^x to be in the region consistent with a -B — >■ r!7 
decay. The selection criteria for -Btag and -Bsig are optimized for each of the r decay modes, 
because the background levels and the background components are mode- dependent. The 
optimization is done so that the figure of merit s/a/s + n is maximized, where s and n are 
the number of signal and background events in the signal enhanced region with remaining 



rameters are calculated in the cms. Here is the nominal B mass 



1^ . and the Pb is 
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energy in the ECL less than 0.2 GeV, which will be described later, calculated with the signal 
branching fraction of 1.79 x 10~^. For leptonic r decays, the dominant background is BB 
pair events correctly tagged by a semileptonic decay. Therefore loose selection criteria are 
chosen to maintain high signal efficiency: 0.5 < P/ < 2.5 GeV/c, —2.1 < cos9b-d*i < 1-3 
for D* mode or —2.6 < cos6b-di < 1-2 for mode, and 0.3 GeV/c < P*^x- -^^^ hadronic 
T decay modes, there is more background from e~^e~ qq continuum and combinatoric 
D^*^^i background. Tighter criteria are used to reduce such backgrounds: 1.0 < P/ < 2.2 
GeV/c, —1.1 < cos ^5_£)(.)0£ < 1.1, and 1.0 < P*^x < 2.4 GeV/c. The upper bound on 
P*_^x is introduced to reject two body B decays. In addition, we suppress continuum back- 
ground by requiring the cosine of the angle between the signal side pion track and the thrust 
axis of the Ptag cos^thr to be less than 0.9. 

The most powerful variable for separating signal and background is the remaining energy 
in the ECL, denoted Pecl, which is the sum of the energies of ECL clusters that are not 
associated with particles from the Ptag and Pgig candidates. We require a minimum energy 
threshold of 50 MeV for the barrel and 100 MeV for the forward, and 150 MeV for the 
backward endcap ECL. A higher threshold is used for the endcap ECL because the effect 
of beam background is more severe. For signal events, Pecl must be either zero or a 
small value arising from residual beam background hits, therefore, signal events peak at low 
Pecl- On the other hand, background events are distributed toward higher Pecl due to the 
contribution from additional particles. We select candidate events in the range Pecl < 1.2 
GeV for further analysis. 

The dominant background contributions are from BB pair and continuum processes, 
which are estimated to be 86% (59%), 11% (36%) for the leptonic (hadronic) r decays from 
MC, respectively. Background from rare P decays and r pair events are estimated from MC 
to be small (3% for leptonic r decays and 4% for hadronic r decay) and fixed in the fit to 
the MC expectation. In order to validate the Pecl simulation, we use a control sample of 
double tagged events, where the Ptag is reconstructed in a semileptonic decay as described 
above and Pgig is reconstructed in the decay chain, P~ —>■ D*^i~u (P*° —> P°7r°), followed 
by — > K~7i~^. The dominant source affecting the Pecl distribution in the control sample 
is beam background hits and is common to the signal. Figure [T] shows the Pecl distribution 
in the control sample for data and the MC simulation scaled to the same luminosity. The 
background in this control sample is negligibly small. The agreement between data and MC 
demonstrates the validity of the Pecl simulation in the signal MC. The Pecl background 
shape is validated using the Pecl sideband region defined by 0.4 < Pecl < 1.2 GeV and 
off-resonance data. We confirm the background shapes obtained from MC agree with the 
Pecl shape of these background data samples. 

After finalizing the signal selection criteria, the signal region is examined. The number of 
signal events is extracted from an extended maximum likelihood fit to the Pecl distribution. 
Probability density functions (PDFs) for each r decay mode are constructed from the MC 
simulation. We use Pecl histograms obtained from MC samples for each of the signal and 
the background components. The PDFs are combined into a likelihood function, 

^=^^nE^./.(^^) (2) 

1=1 j 

where Pj is the Pecl value in the ith event, is the total number of events in the data, 
and Uj is the yield of the jth component, where j is an index for the signal and background 
contributions. In the final fit, four parameters are fioated: the signal yield and the sum of 
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FIG. 1: -Eecl distribution for double semileptonic tagged events. The points with error bars are 
for data and the sohd histogram is the MC expectation. 
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TABLE I: Results of the fit for signal yields and branching fractions. 



BB and continuum backgrounds for the three r decay modes. We combine r decay modes 
by constraining the ratios of the signal yields to the ratio of reconstruction efficiencies 
obtained from MC. Figure [2] shows the E-ecl distribution obtained with the fit results. -Eecl 
distributions for each r decay mode is also shown. We see a clear excess of signal events in 
the region near -Eecl ~ 0. We obtain the signal yield to be ris = I54I35 . The branching 
fraction is calculated as i3 = ns/{2 ■ e ■ Nb+b-), where e is the reconstruction efficiency 
including the tagging efficiency and the branching fraction of the r decay modes and Nb+ b- 
is the number of T(4S') B^B^ events, assuming Nb+b- = ^b^b"' Table [T] shows the 
signal yields and the branching fractions obtained from separate fits for each r decay mode 
and fits with all three modes combined. 

Systematic errors for the measured branching fraction are associated with the uncertain- 
ties in the signal yield, efficiencies and number of B^B^ pairs. The systematic errors for 
the signal yield arise from the uncertainties in the PDF shapes for the signal {Ill'^Vo) and for 
the background (tii;2%) which are dominated by MC statistics. For the latter, uncertainties 
in the branching fractions of B decay modes that peak at -Eecl = such as B~ D'^i'^u 
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with — »• K^TT^^ ^L-^L ^^"i so (-8!4)%! well as uncertainties in the background from 
rare B decays and r pair events (3.8%) are also taken into account. We take a 11.6% error 
as the systematic error associated with the tag reconstruction efficiency from the difference 
of yields between data and MC for the control sample. This value includes the error in the 



branching fraction B{B^ D*^i~i?), which we estimate from B{B^ D*^i^u) in Ref. 14 
and isospin symmetry. The systematic error in the signal efficiencies arises from the un- 
certainty in tracking efficiency (1.0%), particle identification efficiency (1.3%), branching 
fractions of r decays (0.4%), and MC statistics (0.9%). The systematic error due to the 
uncertainty in Nb+b- is 1.4%. The total fractional systematic uncertainty is ^22%) S'^d the 
branching fraction is 

B{B- T-Ur) = (1.65l°J?(stat)l°J5(syst)) x 10"^ (3) 



The significance of the observed signal is evaluated by E = 1/— 2 ln(£o/'^max) where C 
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FIG. 2: Eecl distribution of semileptonic tagged events with the fit result for (a) all r decay modes 
combined, (b) t~ — > e~T'ei^T-, (c) t~ ^~v^i't- and (d) r~ — > 'k~Vt-. The points with error bars 
are data. The hatched histogram and solid open histogram are the background and the signal, 
respectively. 
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and Co denote the maximum likelihood value and likelihood value obtained assuming zero 
signal events, respectively. The systematic uncertainty is convolved in the likelihood with a 
Gaussian distribution with a width corresponding to the systematic error of the signal yield. 
We determine the significance of the signal yield to be 3.8. 

In summary, we have measured the decay t^V with BB pair events tagged by 

semileptonic B decays using a data sample containing 657 x 10^ BB pairs collected at the 
T(4S') resonance with the Belle detector at the KEKB asymmetric e~^e~ collider. We mea- 
sure the branching fraction to be (1.65to;37(stat)lo!37(syst)) x 10"'^. with a significance of 3.8 
standard deviations including systematics. We confirm the evidence reported in the previous 
Belle measurement with BB pair events tagged by hadronic B decays. Using the measured 
branching fraction and known values of G^, rriB, rrir and tb fl^, the product of the B meson 
decay constant Jb and the magnitude of the Cabibbo-Kobayashi-Maskawa matrix element 
\Vub\ is determined to be /b|V^6| = (9.7 ± l.ltij) x 10"^ GeV. The measured branching 
fraction is consistent with the SM expectation from other experimental constraints [5]. 
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